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The weak magnetic field induced reduction of the low-temperature superfluid density in cuprate 
superconductors is studied based on the kinetic energy driven superconducting mechanism. The 
electromagnetic response kernel is evaluated by considering both couplings of the electron charge and 
electron magnetic momentum with a weak magnetic field and employed to calculate the superfluid 
density, then the main features of the weak magnetic field induced reduction of the low-temperature 
superfluid density are qualitatively reproduced. The theory also shows that the striking behavior of 
the weak magnetic field induced reduction of the low-temperature superfluid density is intriguingly 
related to both depairing due to the Pauli spin polarization and nonlocal response in the vicinity of 
the d-wave gap nodes on the Fermi surface to a weak magnetic field. 
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The superfluid density p s , being proportional to the 
density of the supercarriers, is one of the important char- 
acteristic of the superconducting (SC) condensate 2 . It is 
sensitive to the low-lying excitation spectrum, and there- 
fore the knowledge of the superfluid density is essential to 
understanding the physics of the underlying mechanism 
responsible for superconductivity- 2 -^. Since cuprate su- 
perconductors are doped Mott insulators with the strong 
short-range antiferromagnetic correlation dominating the 
entire SC phased, the magnetic field can be also used to 
probe the doping and momentum dependence of the SC 
gap and spin structure of the Cooper pair—. This is why 
the first evidence of the d-wave Cooper pairing state in 
cuprate superconductors was obtained from the experi- 
mental measurement for the magnetic field penetration 
depth A (then the superfluid density p s = A~ 2 ) 5 . Ex- 
perimentally, by virtue of systematic studies using the 
muon-spin-rotation measurement technique, some essen- 
tial features of the superfluid density in cuprate super- 
conductors have been established now for all the tem- 
perature T < T c throughout the SC dome^— . How- 
ever, there are numerous anomalies, which complicate 
the physical properties of the superfluid density. Among 
these anomalies is the magnetic field dependence of the 
superfluid density first observed on the cuprate super- 
conductor YBa 2 Cu 3 07_^— , where a weak magnetic 
field can induce an reduction of the superfluid density 
at the low temperatures. Later, this weak magnetic 
field induced reduction of the low-temperature super- 
fluid density was also found in other families of cuprate 
superconductors^—. 

The appearance of the weak magnetic field induced 
reduction of the low-temperature superfluid density in 
cuprate superconductors is the mostly remarkable effect, 
however, its full understanding is still a challenging issue. 
The earlier work gave the main impetus for a phenomeno- 
logical description of the magnetic field dependent super- 
fluid density in the Meissner state, where it has been ar- 
gued that the weak magnetic field induced reduction of 
the low-temperature superfluid density arises from non- 



linear response of the d-wave state to a weak magnetic 
fields. Later, a more dominant contribution to the mag- 
netic field dependence of the single-particle excitations in 
the superfluid density comes from the nonlocality of the 
supercurrent response in the vicinity of the d-wave gap 
nodes on the Fermi surface in the d-wave SC stated. In 
our recent worki£ based on the kinetic energy driven SC 
mechanism^, the doping and temperature dependence 
of the electromagnetic response in cuprate superconduc- 
tors has been discussed by considering the coupling of 
the electron charge with a weak magnetic field, where 
the Meissner effect is obtained for all the temperature 
T < T c throughout the SC dome, and then the main 
features of the doping and temperature dependence of 
the local magnetic field profile, the magnetic field pen- 
etration depth, and the superfluid density are well re- 
produced. In particular, it is shown that in analogy to 
the domelike shape of the doping dependent SC transi- 
tion temperature, the maximal superfluid density occurs 
around the critical doping 6 w 0.195, and then decreases 
in both lower doped and higher doped regimes. However, 
the coupling of the electron magnetic momentum with a 
weak magnetic field in terms of the Zeeman mechanism 
has been dropped in these discussions 1 -. In this paper, 
we study the weak magnetic field induced reduction of 
the low-temperature superfluid density in cuprate super- 
conductors by considering both couplings of the electron 
charge and electron magnetic momentum with a weak 
magnetic field. Following the linear response theory, we 
have evaluated the magnetic field dependence of the re- 
sponse kernel within the kinetic energy driven SC mech- 
anism. This response kernel is employed to calculate the 
superfluid density, then the main features of the weak 
magnetic field induced reduction of the low-temperature 
superfluid density are qualitatively reproduced. Our re- 
sults also show that the striking behavior of the weak 
magnetic field induced reduction of the low-temperature 
superfluid density is intriguingly related to both depair- 
ing due to the Pauli spin polarization and nonlocal re- 
sponse in the vicinity of the d-wave gap nodes on the 
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Fermi surface to a weak magnetic field. 

We start from the t-J model on a square lattice 2 ^. 
However, for discussions of the weak magnetic field in- 
duced reduction of the low-temperature superfluid den- 
sity in cuprate superconductors, the t-J model can be 
extended by including the exponential Peierls factor and 
Zeeman term as, 

Ifjff 

li) 1 a la 

+ j£SrS^-e*5>c£Ar, (1) 
Ifj la- 
supplemented by an important on-site local constraint 
J2 a C^Cio- < 1 to remove the double occupancy, where 
fj = ±x, ±y, fj' = ±x ± y, C\ a (C/cr) is the electron cre- 
ation (annihilation) operator, S; = (Sf,Sf,Si) are spin 
operators, and fi is the chemical potential. The exponen- 
tial Peierls factors account for the coupling of the elec- 
tron charge to a weak magnetic field in terms of the vec- 
tor potential A(Z), while the Zeeman magnetic energy 
£b = 9HbB accounts for the coupling of the electron 
magnetic momentum g/i^ with the weak magnetic field 
B = rotA, with the Lande factor g and Bohr magne- 
ton jiB- For a proper description of the electron single 
occupancy local constraint, the charge-spin separation 
(CSS) fermion-spin theory 2 ^ has been proposed, where 
the physics of no double occupancy is taken into account 
by representing the electron as a composite object cre- 
ated by Cif — hj^S[~ and C;j. = hj^S^ , with the spin- 

ful fermion operator hi a = e~ l ®'°hi that describes the 
charge degree of freedom of the electron together with 
some effects of spin configuration rearrangements due 
to the presence of the doped hole itself (charge carrier), 
while the spin operator Si represents the spin degree of 
freedom of the electron, then the electron single occu- 
pancy local constraint is satisfied in analytical calcula- 
tions. In this CSS fermion-spin representation, the t-J 
model |T]) can be expressed as, 

H = tY,^ Am (^h ir S+Sr^ + hl^SrS^) 
Ifj 

- *Y, e ~ iiMlw ( h U't h ^sr +fi , +hl n ,^SfS+ +fj ,) 
Ifj' 

la If] I 

where J e g = (1 — S) 2 J, and S = (h, a hia) — {h\h{) is the 
doping concentration. 

For a microscopic description of the SC state of 
cuprate superconductors, the kinetic energy driven SC 
mechanism^ has been developed based on the CSS 
fermion-spin theory, where the charge carrier-spin inter- 
action from the kinetic energy term in the t-J model 



([2]) induces a charge carrier d-wave pairing state by ex- 
changing spin excitations in the higher power of the dop- 
ing concentration, then the electron Cooper pairs orig- 
inating from the charge carrier pairing state are due 
to the charge-spin recombination, and their condensa- 
tion reveals the d-wave SC ground-state. In particular, 
this kinetic energy driven SC state is the conventional 
Bardeen-Cooper-Schrieffer (BCS)-like with the d-wave 
symmetry, so that all main low energy features of the SC 
coherence of the quasiparticle peaks have been quantita- 
tively reproduce d 2 ^ 23 , although the pairing mechanism 
is driven by the kinetic energy by exchanging spin excita- 
tions. For discussions of the weak magnetic field induced 
reduction of the low-temperature superfluid density, we 
generalize the analytical calculation from our previous 
casein without considering the coupling of the electron 
magnetic momentum with a weak magnetic field to the 
case in the presence of the coupling of the electron mag- 
netic momentum with a weak magnetic field. Following 
the discussions in Ref. Qjl and Ref. HH, the mean-field 
(MF) spin excitation spectrum in the t-J model © in 
the presence of the coupling of the electron magnetic mo- 
mentum with a weak magnetic field can be evaluated as 

Wp B ' = y^ojp + (2eb) 2 , with Up is the spin excitation 

spectrum in the case without the coupling of the elec- 
tron magnetic momentum with a weak magnetic field, 
and has been given in Ref. [22T Obviously, an additional 
spin gap 2eb = ^g/isB in the spin excitation spectrum 
is induced by the externally applied magnetic field B. In 
this case, the full charge carrier Green's function can be 
obtained explicitly in the Nambu representation as, 

IT fir i, , R\ - lUnT ° + & T3 ~ ^S^th fos 
Mk, luJ n , B) - Z hF — -j^- , (3) 

K) 2 - Kk 

where To is the unit matrix, t\ and T3 are Pauli ma- 
trices, the renormalized charge carrier excitation spec- 
trum £k = ZhF£k, with the MF charge carrier excita- 
tion spectrum £k = Ztxijk ~ Zt'x2j] i — A*, the spin 
correlation functions \i = {SfSr^.}, X2 = {S+S~ +fr }, 

7k = (l/^E^, 7k = (WEfl'* ^', Z . is the 
number of the nearest neighbor or next-nearest neighbor 
sites, the renormalized charge carrier d-wave pair gap 

_ ( b\ ( B) ~ ( B) 

Aj lZ ; (k) = (k), with the effective charge car- 

rier d-wave pair gap A^ ; (k) = A^ '(cosk x — cosfc w )/2, 

( B) 

and the charge carrier quasiparticle spectrum = 
y^k + |Aj^(k)| 2 , while the magnetic field dependence 
of the effective charge carrier pair gap A^ (k) and 

( B) 

the quasiparticle coherent weight Z^ F satisfy the fol- 
lowing two equations A^(k) = £3 (k, w = 0, B) 
and Zgp-\B) = 1 - E^(k,w = 0,B) | k=M] , 

where (k, w, B) and E^ (k, u, B) are the charge car- 
rier self-energies obtained from the spin bubble, while 
(k,oj,B) is the antisymmetric part of E^ (k, oj,B). 
Moreover, the forms of the obtained E^ (k, u>,B) and 
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Eg (k, (j, B) in the present case are the same as these 
given in Ref. [H for the case without considering the 
coupling of the electron magnetic momentum with a 
weak magnetic field except the spin excitation spectrum 

( B) 

w p has been replaced by u p , therefore the effect of 
a weak magnetic field B on the SC state is consid- 
ered explicitly through the weak magnetic field B enter- 
ing in E^ (k, lj, B) and Eg ^(k, w,B). These two equa- 
tions should be solved simultaneously with other self- 
consistent equation s 19 ' 22 . However, in the limit of the 
weak magnetic field B <^ B c for a given doping concen- 
tration, the main role of the coupling of the electron mag- 
netic momentum with a weak magnetic field is induced 
an additional spin gap in the spin excitation spectrum as 
we have mentioned above, while the effect on other spin 
correlation functions is negligible. In this case, we will 
concentrate on the effect of a weak magnetic field on the 
charge carrier part since the SC state is dominated by 
the charge carrier d-wave pairing state, where the effec- 
tive charge carrier gap parameter ; , the quasiparticle 

( B) 

coherent weight Z^ F , the chemical potential, and other 
charge carrier particle-hole parameters must be solved 
self-consistently. 

In the presence of the coupling of the electron magnetic 
momentum with a weak magnetic field, the magnetic field 
dependence of the response current density J ^ and the 
vector potential A v are related by a nonlocal kernel of 
the response function K^ v asi£, 

3 

J M (q,a;,B) = - K ^(^ w, B)A v (q, u), (4) 
i/=i 

with the Greek indices label the axes of the Cartesian 
coordinate system. This magnetic field dependence of 
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4 j (k,q) 



It is easy to show^ that in the long wavelength limit, 
i.e., |q| — > 0, K$ (q — > 0,0, B) = at the temperature 
T = 0. In this case, the long wavelength electromagnetic 
response is determined by the diamagnetic part of the 
kernel only. On the other hand, at the SC transition tem- 
perature T = T c , K$(q -> 0,0,5) = -(1/A|), which 
exactly cancels the diamagnetic part of the response ker- 
nel (pjaj) . and then the Meissner effect in the presence of 
the coupling of the electron magnetic momentum with a 
weak magnetic field is obtained for all T <T C throughout 



the response kernel ([4| can be separated into two parts 
as K^ v (q,u,B) = K$(q,u>,B) + K$ (q, uj , B) , where 
Kffi(q,uj,B) and K$ (q,w, B) are the corresponding 
diamagnetic and paramagnetic parts, respectively. In the 
CSS fermion-spin representation 2 ^, the vector potential 
A has been coupled to the electron charge, which are 
now represented by C/-f = hj^.Sf and C;j, — hj,S~^. In 
this case, the electron polarization operator is expressed 
as P = — eJ2 R-iCiijCia — e^2~R,ih\hi, then the corre- 

■ia ia 

sponding electron current operator is obtained by eval- 
uating the time-derivative of this polarization operator 
as j = + j j with and j( p ) are the correspond- 
ing diamagnetic (d) and paramagnetic (p) components 
of the electron current operator. Following our previous 
discussions^, these diamagnetic and paramagnetic parts 
of the magnetic field dependence of the response kernel 
KjfJ(q,uj,B) and i4T^(q, w, B) can be obtained in the 
the static limit as, 

K$(q,0,B) = -jLfafat -2x2^)8^ = i^, (5a) 
K$(q,0,B) = l^ 7A1 (k + q,k) 7 :(k + q,k)[4 S) (k,q) 

k 

+ L 2 B \k,q)]=K^(q,0,B)S^, (5b) 

where the charge carrier particle-hole parameters 4>\ — 
(hj a h l+fia ) and 2 = (/jt^i-Myv), A^ 2 = -Ae 2 {xi4>it - 
2x2<f>2t')/fi 2 is the London penetration depth, and now 
is doping, temperature, and magnetic field dependent, 
the bare current vertex 7 M (k + q, k) has been given in 

Ref. [H, while the functions L[ B \k., q) and Zi£'(k, q) 
are obtained as, 



(6a) 



(6b) 

I 

the SC dome. 

However, the result we have obtained the response ker- 
nel in Eqs. (|5a| and (j5b[) can not be used for a direct 
comparison with the corresponding experimental data of 
cuprate superconductors because the response kernel de- 
rived within the linear response theory describes the re- 
sponse of an infinite system, whereas in the problem of 
the penetration of the field and the system has a surface, 
i.e., it occupies a half-space x > 0. In such problems, 
it is necessary to impose boundary conditions for charge 
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carriers. This can be done within the simplest specular 
reflection model 2 ^ with a two-dimensional (2D) geometry 
of the SC plane. Taking into account the 2D geometry 
of cuprate superconductors within the specular reflection 
model 2 ^, wei 8 - can obtain the magnetic field penetration 
depth as, 



A(T, B) = i J h z (x,B)dx 
o 

oo 

dq x 



7T J jJioKyy(q x ,Q,Q,B) + ql' 



(7) 



which therefore reflects the measurably electromagnetic 
response in cuprate superconductors. 
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FIG. 1: The magnetic field penetration depth AA as a func- 
tion of temperature at the doping concentration 8 = 0.09 for 
the magnetic field B = T (solid line), B — 0.5 T (dashed 
line), and B = 1.0 T (dash-dotted line) with parameters 
t/J = 2.5, t'/t = 0.3, and J = 1000K. Inset: the corre- 
sponding experimental data for YBa2Cu4 0s taken from Ref. 

Now we are ready to discuss the magnetic field depen- 
dence of the Meissner effect in cuprate superconductors. 
In cuprate superconductors, although the values of J, t, 
and t' are believed to vary somewhat from compound 
to compound^, however, as a qualitative discussion as in 
our recent worki 8 -, the commonly used parameter a 25 i 26 
in this paper are chosen as t/J = 2.5, t'/t = 0.3, 
and J = I000K. Furthermore, a characteristic length 
scale ao = yj h 2 a/ /ioe 2 J is introduced. Using the lat- 
tice parameter a rs 0.383nm for YBa 2 Cu307_ a , this 
characteristic length is obtain as ao ~ 97.8nm. In this 
case, the magnetic field penetration depth AA(T, B) = 
A(T, B) — A(0, 0) as a function of temperature at the 
doping concentration 6 — 0.09 for the magnetic field 
B = T (solid line), B = 0.5 T (dashed line), and 
B = 1.0 T (dash-dotted line) is plotted in Fig. [1] m com- 
parison with the corresponding experimental results^ of 
YBa2Cu40g (inset). The similar magnetic field depen- 
dence of the magnetic field penetration depth has been 
also observed experimentally on YE^CuaOg.gs 8 . How- 
ever, YBa2Cu3 07_,5 contains a single CuO chain per unit 



cell whose oxygen content can be varied to tune the dop- 
ing level on the Cu0 2 planed, while its close relative 
YBa2Cu408 has a double chain layer that is stoichiomet- 
ric and a planar state that is underdopecU^. Although 
this difference of the electronic structure of the quasi- 
2D plane leads to some subtly different behaviors be- 
tween YBa2Cu307_,5 and YBa2Cu40s, where at a weak 
magnetic field, the experimental curves of the temper- 
ature dependent magnetic field penetration depth for 
YBa 2 Cu3 06.95 show curvature in low temperature 8 - as in 
the case of zero magnetic field, while the corresponding 
experimental curves of the temperature dependent mag- 
netic field penetration depth for YBa2Cu40s are linear—, 
the qualitative properties in both YBa2Cu307_a and 
YBa2Cu4C>8 are consistent each other— In this pa- 
per, we mainly focus on the qualitative properties of the 
magnetic field dependent Meissner effect in cuprate su- 
perconductors based on the simple t-J model (JI]). Within 
the kinetic energy driven SC mechanism, the SC transi- 
tion temperature T c = 52K at 8 = 0.09 for zero mag- 
netic field. As wei 8 - have shown that at the SC tran- 
sition temperature T — T c , the kernel of the response 
function K )lv (q — > 0, 0, 0)\t=t c = 0. In this case, we ob- 
tain the magnetic field penetration depth from Eq. ([7]) as 
A(T C , 0) = oo, which reflects that in the normal state, the 
external magnetic field can penetrate through the main 
body of the system, therefore there is no the Meissner 
effect in the normal state. Furthermore, our present re- 
sult in Fig. [T] shows clearly that the characteristic feature 
of the temperature dependent A(T, B) is essentially inde- 
pendent on a weak magnetic field, in particular, A(T, B) 
shows a crossover from the linear temperature depen- 
dence at higher temperatures to a nonlinear one in the 
low temperatures as in the case of zero magnetic field 18 , 
in qualitative agreement with the corresponding exper- 
imental data 8 - of YBa2Cu3 06.95- Moreover, the magni- 
tude of X(T, B) at the low temperatures is dependent on a 
weak magnetic field, and then it is independent on a weak 
magnetic field at higher temperatures. However, there is 
a substantial difference between theory and experiment, 
namely, the value of the magnetic field dependent pen- 
etration depth AA(T, B) at the low temperatures cal- 
culated theoretically is smaller than the corresponding 
value measured in the experiment 8 -^. However, upon 
a closer examination one can see immediately that the 
main difference is due to fact that the calculated value 
of AA(T, B) increases slowly with magnetic field at the 
low temperatures. The simple t-J model can not be re- 
garded as a complete model for the quantitative com- 
parison with cuprate superconductors, however, as for 
a qualitative discussion in this paper, the overall shape 
seen in the theoretical result is qualitatively consistent 
with that observed in the experiment&i^. 

To show this magnetic field dependence of A(T, B) at 
the low temperatures clearly, we have made a series of cal- 
culations for X(T, B) at differently weak magnetic fields, 
and the result of AA(T, B) as a function of magnetic field 
at S = 0.09 with T = 4K is plotted in Fig. OH in com- 
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FIG. 2: The magnetic field penetration depth as a function 
of magnetic field at the doping concentration S = 0.09 for 
temperature T — 4K with parameters t/J — 2.5, t'/t = 0.3, 
and J = 1000K. Inset: the corresponding experimental data 
for YBa2Cu30e.95 at temperature T — 4.2K (circles) and T = 
7K (squares) taken from Ref. [nj. 

parison with the corresponding experimental results^ of 
YBa 2 Cu3 06.95 at temperature T = 4.2K (circles) and 
T = 7K (squares) (inset). Obviously, X(T,B) is a nonlin- 
ear function of magnetic field. Moreover, a rapid increase 
of A(T, B) with a weak magnetic field at the low tem- 
peratures observed from cuprate superconductor s 10 ' 11 is 
qualitatively reproduced. It should be emphasized that 
the present result for the d-wave SC state in cuprate su- 
perconductors is much different from that in the conven- 
tional superconductors, where the magnetic field depen- 
dence is typically weak at the low temperatures because 
the isotropic energy gap exponentially cut off the quasi- 
particle excitations. 
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FIG. 3: The superfluid density as a function of temperature 
at the doping concentration 5 — 0.09 for the magnetic field 
B = 0T (solid line), B = 0.5 T (dashed line), and B = 1.0 T 
(dash-dotted line) with parameters t/J — 2.5, t' /t = 0.3, and 
J = 1000K. Insets: the corresponding experimental results 
for YBa2Cu3 06.95 taken from Ref. [Tct 

Next we discuss the weak magnetic field induced reduc- 



tion of the low-temperature superfluid density p s (T,B). 
This superfluid density is closely related to the magnetic 
field penetration depth as p s (T,B) = \~ 2 (T,B). In this 
case, we have made a series of calculations for p s (T,B) 
at differently weak magnetic fields, and the result of 
p s (T, B) as a function of temperature at S = 0.09 for B = 
T (solid line), B = 0.5 T (dashed line), and B = 1.0 T 
(dash-dotted line) is plotted in Fig. [3]in comparison with 
the corresponding experimental data for YBa2Cu30g.95 
taken from Ref. US1 (inset). In correspondence with the 
result of A(T, B) shown in Fig. [TJ the characteristic fea- 
ture of the temperature dependent superfluid density is 
also independent on a weak magnetic field, where as in 
the case of zero magnetic field, p s (T,B) shows a lin- 
ear temperature dependence at higher temperatures, and 
then it crosses over to a nonlinear temperature behavior 
at the low temperatures. However, most importantly, the 
magnitude of p s {T, B) at the low temperatures decreases 
with increasing magnetic field, and then it turns to be 
independent on a weak magnetic field at higher temper- 
atures, in qualitative agreement with experimental data 
of cuprate superconducto ra 9 i 10 ' 12 ~ — . The present result 
also indicates that the nature of the quasiparticle exci- 
tations at the low temperatures is strongly influenced 
by a weak magnetic field. This weak magnetic field in- 
duced reduction of the low-temperature superfluid den- 
sity in cuprate superconductors contrasts with the that 
observed from the conventional superconductors^!, where 
the curves of the temperature dependent superfluid den- 
sity for differently weak magnetic fields were found to 
collapse onto a single curve since the conventional super- 
conductors are fully gaped. 

The weak magnetic field induced reduction of the low- 
temperature superfluid density in cuprate superconduc- 
tors arises from both depairing due to the Pauli spin po- 
larization and nonlocal response in the vicinity of the 
d-wave gap nodes on the Fermi surface to a weak mag- 
netic field. In the framework of the kinetic energy driven 
SC mechanism, the d-wave SC state is mediated by the 
interaction of electrons and spin excitations^ 9 -, where the 
depairing can occur due to the Pauli spin polarization 
in the presence of an external magnetic field. This fol- 
lows a fact that an applied magnetic field aligns the 
spins of the unpaired electrons, i.e., there is a tendency 
to induce the magnetic order, then the kinetic energy 
driven d-wave Cooper pairs can not take advantage of 
the lower energy offered by a spin-polarized stat o 19 ' 28 . 
On the other hand, the characteristic feature of the 
d-wave superconductors is the existence of four nodes 
on the Fermi surface, where the energy gap vanishes 
A h (k) | at nodes = A h (cosfc 2; -cosfc y )/2| a t nodes = 0. These 
two special features in the kinetic energy driven d-wave 
SC state indicate that even small thermal energy or exter- 
nally small magnetic energy can excite excitations, then 
the superfluid density decreases with increasing tempera- 
ture or increasing magnetic field, reflecting that the weak 
magnetic field induced reduction of the low-temperature 
superfluid density is a natural consequence of the kinetic 
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energy driven d-wave SC state. Since the quasiparticles 
selectively populate the nodal region at the low temper- 
atures, then the most physical properties in the SC state 
are controlled by the quasiparticle excitations around the 
nodes. In this case, the Ginzburg-Landau ratio around 
the nodal region is no longer large enough for the sys- 
tem to belong to the class of type-II superconductors, 
and the condition of the local limit is not fulfilled 1 - ! 18 ! 29 . 
On contrary, the system falls into the extreme nonlo- 
cal limit, then the nonlinear characteristic in the tem- 
perature dependence of the superfluid density (then the 
magnetic field penetration depth) can be observed exper- 
imentally in cuprate superconductors at the low temper- 
atures. However, when a weak magnetic field is applied 
to the system even at zero temperature, the quasiparti- 
cles around the nodal region become excited out of the 
condensate, and at the same time the electron attractive 
interaction for the Cooper pairs by exchanging spin exci- 
tations is weaken^, both these effects lead to a decreases 
in the superfluid density. With increasing temperatures, 
the externally small magnetic energy due to the presence 
of a weak magnetic field is comparable with the small 
thermal energy at the low temperatures, therefore both 
small thermal energy and weak magnetic field induce an 
reduction of the superfluid density. However, at higher 
temperatures, this externally small magnetic energy is 
much smaller than the thermal energy, then the major 
contribution to a decrease of the superfluid density comes 
from the thermal energy. This is why a weak magnetic 
field only reduces an reduction of the superfluid density 
only at the low temperatures. 

As we have mentioned in Eq. (pj , the coupling of the 
electron charge to a weak magnetic field is in terms of the 
vector potential A(7), while the coupling of the electron 
magnetic momentum g/iB with the weak magnetic field 
B = rotA is in terms of the Zeeman mechanism. Dur- 
ing the above discussions, the electromagnetic response 
kernel (O is calculated with the bare current verte»i^, 
where the important point is that for a weak magnetic 
field which orientation is the same at all spatial points 
the spin polarization axis along the weak magnetic field 
is chosen accordingly. As a consequence, we do not take 
into account longitudinal excitations properly^, and the 



obtained results are valid only in the gauge, where the 
vector potential is purely transverse, e.g. in the Coulomb 
gauge. If we want to keep the theory gauge invariant, it is 
crucial to approximate the electromagnetic response ker- 
nel in a way maintaining local charge conservation^. For 
the case without considering the coupling of the electron 
magnetic momentum with a weak magnetic field, we have 
shown that although the electromagnetic response kernel 
is not manifestly gauge invariant within the bare current 
verte»i^, the gauge invariance is kept within the dressed 
current verte?*^. However, for the present case with con- 
sidering both couplings of the electron charge and elec- 
tron magnetic momentum with a weak magnetic field, 
the gauge invariance should be kept within the dressed 
current vertex together with accordingly rotated spin po- 
larization axis along the weak magnetic field. These and 
the related issues are under investigation now. 

In conclusion, we have discussed the magnetic field in- 
duced reduction of the low-temperature superfluid den- 
sity in cuprate superconductors based on the kinetic en- 
ergy driven SC mechanism by considering both couplings 
of the electron charge and electron magnetic momentum 
with a weak magnetic field. Our results show that al- 
though the characteristic feature of the temperature de- 
pendent superfluid density is found to be independent on 
a weak magnetic field, this weak magnetic field induces 
an reduction of the low-temperature superfluid density 
in the Meissner state. Our results also show that the 
striking behavior of the weak magnetic field induced re- 
duction of the low-temperature superfluid density can be 
attributed to both depairing due to the Pauli spin polar- 
ization and nonlocal response in the vicinity of the d-wave 
gap nodes on the Fermi surface to a weak magnetic field. 
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